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Abstract: In Hebei Province of North China, forest was recovered with natural recruitment in plantations 
with large area of clear-cutting Chinese pine (Pinus tabuliformis). This study was aimed to demonstrate the 
dynamic characteristics of recruits during the natural recruitment. Both plot survey and the spatial 
point-pattern analysis were performed. Five developmental stages of natural recruitment were selected 
and studied, including 1 year before and 2, 5, 8, and 11 years after clear-cutting. Different slope aspects 
were also included. Natural recruitment was always dominated by Chinese pine with a proportion of 
higher than 90%. For plots of 1 year before clear-cutting on east- and north-facing slopes, recruit densities 
were 7886 and 5036 stems/hm?, the average heights were 0.78 (£0.85) and 1.06 (+1.15) m, and the 
average diameters at breast height (DBH) were 3.21 (+1.38) and 2.91 (+1.38) cm, respectively. After 
clear-cutting, recruit density was initially increased, then it was gradually declined with time; however, the 
variation of average DBH was contrary to that of recruit density. Both of them were no longer varied 
between 8 and 11 years after clear-cutting. The average height of recruits continued to increase after 
clear-cutting. For the plots of 11 years after clear-cutting on east- and north-facing slopes, average heights 
of recruits reached 2.00 (£1.14) and 2.24 (+1.20) m, respectively. The statuses of recruits on north-facing 
slopes wete better than those on east-facing slopes after clear-cutting. Meanwhile, recruits on east-facing 


slopes were always aggregated at small scales, while spatial pattern of recruits varied with time on 
north-facing slopes. Moreover, forest was recovered more quickly by natural recruitment than by artificial 
afforestation after clear-cutting. The structural diversity was higher in naturally regenerated forests than in 
plantations of the same age. Our results demonstrated that clear-cutting of Chinese pine plantations 
recovered by natural recruitment has the potential to be an effective approach for establishing 
multifunctional forest. 
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1 Introduction 


Chinese pine (Pinus tabuliformis) is an important vegetation species in North China (Wu, 1980). 
From 1950$, several plantations of Chinese pine have been built by government in North China 
for wood production. Now, a large area of Chinese pine plantations has distributed over North 
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China. Since the general felling age of Chinese pine plantations was 50 years, many plantations 
were in the period of final felling. Clear-cutting was generally applied by the management 
departments. Then the forest was expected to be reconstructed with natural recruitment, because 
of many potential advantages, including higher recruit density (Prévost and Gauthier, 2013), 
structural diversity (Yamagawa et al., 2008), better conservation of local vegetation genotypes 
(Peterken, 1996) and reduction of management costs (CCFM, 2006). 

Many studies found that forest could be reconstructed, and the reconstructed forest relying on 
natural recruitment was close to natural forest in conifer plantations after clear-cutting in 
temperate regions (Yamagawa and Ito, 2006; Man et al., 2013; Spracklen et al., 2013; Fang et al., 
2014). Thus, the further investigation on the process of natural recruitment could be conductive to 
predicting succession tendency and optimizing management strategies (Sakai et al., 2006; Swaim 
et al., 2016). Until 2015, the earliest clear-cutting Chinese pine plantations have been restored for 
11 years in North China. However, as we know, no related literatures have been reported on the 
recruitment dynamics of Chinese pine. 

In general, current studies on the recruitment after clear-cutting have mainly performed from 
the following aspects: composition of tree species (Coates, 2002; Haeussler and Bergeron, 2004; 
Taylor et al., 2013), vertical structure of recruits (Buermeyer and Harrington, 2002; Wu et al., 
2008; Spracklen et al., 2013), and characteristics of the recruitment, including density (Sakai et al., 
2006; Ilisson and Chen, 2009; Man et al., 2013), growth status (Stuiver et al., 2016; Swaim et al., 
2016), and age (Sturgess and Atkinson, 1993; Buermeyer and Harrington, 2002; Kurokochi et al., 
2010). Moreover, researchers also found that the dynamics of natural recruitment was related to 
the spatial pattern of vegetation among different heights and age classes (Fajardo et al., 2006; 
Franklin and Santos, 2011; Luo et al., 2012; Carrer et al., 2013). Although all these aspects were 
relevant to the recruitment after clear-cutting, two critical factors were considered in our study for 
further investigation. On the one hand, in the process of forest restoration after clear-cutting, 
advanced regeneration was proved to be an essential part of the forest land vegetation (Buermeyer 
and Harrington, 2002; Greene et al., 2002; Yamagawa and Ito, 2006). In addition, the older 
advanced recruit was an important seed source for further forest reconstruction (Ilisson and Chen, 
2009). Thus, we paid attentions to the advanced regeneration of 1 year before clear-cutting. It was 
aimed to reveal the growth conditions of advanced regeneration, to explore the effects of 
clear-cutting on advanced regeneration, and to find the relations between advance regeneration 
and forest restoration. On the other hand, slope aspect was one of the dominating factors 
influencing the natural recruitment of Chinese pine plantations (Xu, 1993; Ren et al., 2010; Han 
et al., 2012; Liu et al., 2014). Therefore, the effects of slope aspect on characteristics of natural 
recruitment were also analyzed in this study. As mentioned above, the aims of the study were (1) 
to elucidate the dynamic of advanced regeneration and its role on the process of forest recovery 
after clear-cutting; (2) the variation trends of natural recruitment over time after clear-cutting; and 
(3) the effects of slope aspect on natural recruitment. 


2 Materials and methods 


2.1 Study area 


The study area of Chinese pine plantations was situated in Pingquan County, Hebei Province, 
North China (Fig. 1). It was approximately 250 km northeast of Beijing. It belongs to the 
temperate continental monsoon mountain climate. The annual mean temperature was 7.4°C and 
the annual precipitation was 356.7—620.3 mm from 2004 to 2015 (Table 1). The main soil type 
was mountain cinnamon soil with gravel. In the study area, the area of plantations has accounted 
for 95% of the total forest area in Pingquan County, with four main species including Chinese 
pine, aspen (Populus davidiana), Mongolian oak (Quercus mongolica) and black locust (Robinia 
pseudoacacia). And the area of Chinese pine plantations was the largest, with an area of 8901.85 
hm? and accounted for 61.6% of the local total forest area. 
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Fig.1 Sampling plots and study area 


Table 1 Seasonal and annual precipitation (mm) of the study area 


Year Spring Summer Autumn Winter Full-year 
2004 118.6 59.5 316.5 15.9 510.5 
2005 27.4 106.2 444.4 9.1 587.1 
2006 47.2 36.8 330.3 SM) 420.2 
2007 85.6 140.6 303.1 6.8 536.1 
2008 114.3 90.9 307.7 3.9 516.8 
2009 101.3 68.0 167.3 20.1 356.7 
2010 203.1 86.0 322.0 92 620.3 
2011 70.6 552 388.8 8.2 522.8 
2012 100.4 70.3 334.1 9.5 514.3 
2013 81.6 37.6 354.3 4.0 477.5 
2014 90.0 91.3 280.4 25 464.0 
2015 128.2 98.6 297.1 8.5 532.4 


Note: Spring means the period of March-May; Summer, June-August; Autumn, September-November; Winter, December—February of 
the next year. 


2.2 Data collection 


The samples were taken along the chronosequence. Ten forest stands were selected with different 
developmental stages and slope aspects. Two slope aspects were involved: east and north (л=5 for 
each). Chinese pine plantations were mainly distributed on half shady (east- and northwest-facing) 
and shady (north- and northeast-facing) slopes due to the limitation of precipitation in this region 
(Table 1). Thus, in this study, east- and north-facing slopes were selected to represent half shady and 
shady slopes, respectively. Five forest stands on east-facing slopes were taken from plantations of 1 
year before clear-cutting, and 2, 5, 8, and 11 years after clear-cutting, which were named as EB1, 
EA2, EAS, EAS, and EA11, respectively. The other 5 forest stands on north-facing slopes were also 
taken from plantations of 1 year before clear-cutting, and 2, 5, 8 and 11 years after clear-cutting, 
which were named as NB1, NA2, NA5, NA8, and NATI, respectively. The characteristics of forest 
stands were shown in Table 2. For all the plantations, Chinese pine was the only one canopy wood 
species and there was also only one canopy layer. The stands EB1, NB1 were 49-years old Chinese 
pine plantations. However, the general cutting age of Chinese pine plantations was 50, thus the 
stands EB1, NBI were applied to represent the Chinese pine plantations 1 year before clear-cutting. 
The sampling of stands was aimed to represent the status of advanced regeneration before 
clear-cutting and demonstrated the effects of clear-cutting on advanced regeneration. For the 
plantations 2, 5, 8 and 11 years after clear-cutting, all canopy woods were removed by clear-cutting. 
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Meanwhile, no artificial afforestation was carried out after clear-cutting. The forest recovery was 
completely depended on natural recruitment. 

The growth characteristics of vegetation lasted for years and were hard to be continuously tracked. 
A few studies reported the variation trends of vegetation after clear-cutting (Bermúdez et al., 2007). 
Therefore, comparing the growth characteristics of vegetation in different growth stages was a 
convenient method for exploring long-term variations (Bermudez et al., 2007; Akhavan et al., 2012; 
Stuiver et al., 2016). 

The field sampling and investigation were performed from June to September in 2015. The 
selected plots were in the center of each forest stand, with an area of 60 mx60 m. The distance from 
plot boundary to edge of stand or cutting blank was at least 50 m. Only one plot was sampled at 
each stand, which was named after corresponding forest stands. For each plot, the species, quantity, 
height, coordinates of all recruits (with the height 20.1 m) and the DBH of recruits (with the height 
22 m) were recorded. For the plots in stands EB1 and NBI, the coordinates of canopy trees were 
also recorded. In addition, ten quadrates with the area of 5 mx5 m were randomly set in every plot. 
For each quadrate, the age of all recruits was also recorded. According to height (h), we divided the 
recruits into seedlings (0.1 п1</<0.5 m) and saplings (h>0.5 m). Only seedlings higher than or equal 
to 0.1 m were included because only the established recruitment was considered (i.e., seedlings that 
survived at least 2 years) (Ignacio et al., 2009). Most of the observed recruits were Chinese pine, 
with high mortality rate in the first two years (Xu, 1993). The age of recruits was analyzed on a 
tree-ring analysis system with increment cores or basal sections (WinDENDROTM, 2005). 


2.3 Data analysis 


2.3.1 ANOVA test 


One-way ANOVA test was applied to analyze the differences of height, DBH, and age of recruits 
among plots with different slope aspects or developmental stages. The data analysis was performed 
following a randomized complete block design. The normality and equal variance of data were 
checked and an unequal variance model was applied if necessary. 


2.3.2 Spatial point pattern analysis 


In the spatial point pattern analysis, all information was included in Ripley's K function with a 
certain radius (scale). However, with the increased radius (scale), the information on small distance 
would be included in the results on long distance (scale). The effects of small and large scales would 
be confused with this cumulative calculation (Condit et al., 2000; Schurr et al., 2004; Wiegand and 
Moloney, 2004). 

Moloney's O-ring statistical methods were based on Ripley's K function and Mark correlation 
function. However, the ring of width w was applied in Moloney's O-ring statistical method instead 
of the circle of radius r in Ripley's К function. The average ring point number in the ring was 
calculated. Thus, the special distance grade was isolated. Therefore, Ripley's K function can be 
applied to distinguish distribution pattern of study object beyond a specific scale r. O-ring statistical 
method was based on probability density function, so the distribution pattern on specific scale r 
could be more intuitively reflected (Wiegand and Moloney, 2004). Both the univariate statistics and 
bivariate statistics were included in O-ring statistics. Univariate statistics was applied to analyze the 
distribution pattern of a single object represented by O(r) statistics, while, the bivariate statistics 
was applied to analyze the distribution pattern of the two objects represented by O12(r). Statistics of 
univariate and bivariate O-ring were calculated as follows: 


LV Points [Re] 1 Y points IRO] 
o(r)- 5 0005 Os (r= ZO 
I V Area[ R7 (r)] — У Areal Ry; (7)] 


1 i=l Nn j= 


In Equation 1, nı is the number of individuals for one population in the study area, АҮ (т) is a 
ring around the i^ point of this population, with radius of r and width of w, Points)[Ri(r)] is the 
number of individuals for this population in ring, Агеа[А" (7)] is the area of ring (m°). In Equation 
2, nı is the number of individuals for population 1 in the study area, № (г) is a ring around the i‘ 
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point of the population 1, with radius of r and width of w, Points:[R7;(r)] is the number of 
individuals for population 2 in ring and Area[R?(r)] is the area of ring (mô. 

In this study, overall recruits (including both seedlings and saplings) were analyzed with 
univariate O-ring statistical method. The spatial correlations between recruits and canopy woods as 
well as seedlings and saplings were analyzed with bivariate O-ring statistical method. The O-ring 
statistical method required careful selection of null hypothesis (Wiegand and Moloney, 2004; López, 
2010). For univariate analysis of recruits, complete spatial randomness (CSR) was applied as null 
model which was the simplest and most widely used null model for this kind of analysis (Wiegand 
and Moloney, 2004). This null model considered that any point of the pattern has an equal 
probability of occurring at any position in the plot, and the position of a point is independent of the 
position of any other point (i.e., points do not interact with each other) (Wiegand and Moloney, 
2004). 

For bivariate analysis, there were two kinds of situations. When different age classes (recruits and 
canopy woods) were spatially analyzed, we effectively assumed that both age classes represented 
sequential events in time, therefore the development of the pattern of the canopy woods would not 
be influenced by the pattern of recruits, but the opposite may be true (i.e., the development of the 
recruits may be influenced by the canopy woods) (Fajardo and González, 2009). In this case, the 
choice of an appropriate null model was influenced by antecedent conditions (Wiegand and 
Moloney, 2004). For the situation in our study, it meant that for investigating the relationship 
between recruits and canopy woods, an appropriate null model for testing repulsion or attraction 
would be the randomization of recruit locations and kept the canopy woods locations fixed. The 
reason was that the canopy woods were already presented and without changing their positions 
during the establishment of recruitment and the recruits could potentially be found at the entire plot 
(Wiegand and Moloney, 2004; Carrer et al., 2013). Thus, antecedent condition null model was 
applied to test whether formation process of pattern of recruits was influenced by pattern of canopy 
woods. When seedlings and saplings were compared, independence null model was applied. This 
null model assumed that the spatial patterns of seedlings and saplings were generated by two 
independent and stochastic processes (e.g., different dispersion process periods) (Wiegand and 
Moloney, 2004; Fajardo et al., 2006). 

These analyses were performed with software Programita. The statistical significance of 
departures of O(r) from randomness was tested with Monte Carlo simulations and following three 
null models were applied: CSR, antecedent condition, and independence. The confidence limits 
(envelopes) were obtained from the highest and lowest values of the functions, which were taken 
from 99 simulations of the null models. 


3 Results 


3.1 Composition of tree species 


In all plots, the tree species were analyzed and the dominant species was Chinese pine (Table 3). 
The proportions of Chinese pine were greater than 90% in all plots except for NB1. The 
composition of accompanying tree species was different before and after clear-cutting (Table 3). 
Shade-tolerance species Mongolian oak and hawthorn (Crataegus pinnatifida) only existed in 
ЕВ1 and NBI plots. After clear-cutting, more shade-intolerance pioneer species aspen was 
observed. Slope aspects made little effect on the composition of tree species (Table 3). 


3.2 Density of recruits 


The densities of advanced recruits were 7886 and 5036 stems/hm? for ЕВ1 and NBI (1 year 
before clear-cutting), respectively. For the plots after clear-cutting, the density of recruits was 
significantly increased at the first 2 years, then gradually decreased and no longer varied between 
8 and 11 years after clear-cutting (Fig. 2). Density of recruits was significantly different between 
east- and north-facing slopes (Fig. 2). The recruit densities in plots on east-facing slopes were 
higher than those on north-facing slopes (Fig. 2). The recruit density was increased by 
15.90%-56.59%. Specially, the recruit density in NA5 plot was 18.04% higher than that in EAS 
plot (Fig. 2). 
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Table3 Composition of tree species in different plots (stems/plot) 


Treatment Plot Pt От Ср As Pd Up Rp Tt Bp Sm 
1 year before EB1 2725 57 16 21 1 17 2 

clear-cutting МВ] 1621 47 74 27 12 32 

2 years after EA2 5259 149 14 

clear-cutting NA2 4299 4 54 8 4 
5 years after EA5 2533 24 71 12 

clear-cutting NAS 3045 1 145 17 13 

8 years after EA8 1810 23 36 

clear-cutting NA8 1488 1 111 9 

11 years after _ БА11 1771 4 Th 10 п 

clear-cutting МА11 1536 2 71 7 


Note: Pt, Pinus tabulaeformis; Qm, Quercus mongolica; Cp, Crataegus pinnatifida; As, Armeniaca sibirica; Pd, Populus davidiana; 


Up, Ulmus pumila; Rp, Robinia pseudoacacia; Tt, Tilia tuan; Bp, Betula platyphylla; Sm, Salix matsudana. 
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Fig. 2 Variation of density of recruits in different plots. B1, 1 year before clear-cutting; A2, 2 years after 
clear-cutting; A5, 5 years after clear-cutting; A8, 8 years after clear-cutting; A11, 11 years after clear-cutting. 


3.3 Height of recruits 


In the plots of 1 year before clear-cutting, the average heights of recruits were 0.78(+0.85) and 
1.06(+1.15) m for ЕВІ and NBI, respectively. After clear-cutting, the average height of recruits 
was gradually increased with time (Fig. 3). The final average heights of recruits were 2.00(+1.14) 
and 2.24(+1.20) m in EA11 and NA11 plots, respectively, which were significantly increased by 
156.0096 and 110.8096 (P«0.01). The average height of recruits would be observably influenced 
by slope aspects (Fig. 3). The average heights of recruits in plots on north-facing slopes were 
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Fig. 3 Variation of height of recruits in different slopes and plots. B1, 1 year before clear-cutting; A2, 2 years 
after clear-cutting; A5, 5 years after clear-cutting; A8, 8 years after clear-cutting; All, 11 years after clear-cutting. 
Bars mean standard errors. 
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significantly higher than those on east-facing slopes (P«0.01 for all; Fig. 3). The distribution of 
recruit height would be significantly influenced by the slope aspect, i.e., the distribution was more 
decentralized on the north-facing slopes (Fig. 4). 


East-facing slope North-facing slope Chronosequence 
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Fig. 4 Distribution of height of recruits in different plots. Height classes: I (0.1—0.5 m), II (0.5—1.0 m), III 
(1.0-1.5 m), IV (1.5-2.0 m), V (2.0-2.5 m), VI (2.5-3.0 m), VII (3.0-3.5 m), VIII (3.5-4.0 m), ІХ (4.0-4.5 m), X 
(4.5-5.0 m), and XI (55.0 m). 


3.4 Diameter at breast height (DBH) of recruits 


In the plots of 1 year before clear-cutting, the average DBH of recruits were 3.21 (+1.38) and 2.91 
(+1.38) cm for EB1 and NBI, respectively. After clear-cutting, the average DBH of recruits was 
firstly decreased, then gradually increased and no longer varied between 8 and 11 years (Fig. 5). 
The average DBH of recruits would be significantly influenced by slope aspects (Fig. 5). Before 
clear-cutting, the average DBH of recruits in EB1 plot was higher than that in NBI plot (P=0.015; 
Fig. 5). However, after clear-cutting, the average DBH of recruits in plots on north-facing slopes 
was higher than that on east-facing slopes (P«0.01 for all; Fig. 5). 


ChinaXiv& f'ERBTHI 


DONG Bogian et al.: Forest recovery after clear-cutting in Chinese pine (Pinus tabuliformis)... 


6 
С] East-facing slope 
E North-facing slope 


0 


Bl A2 А8 А11 
Chronosequence 


o 


DBH (cm) 


N 


Fig.5 Diameter at breast height (DBH) of recruits in different slopes and plots. B1, 1 year before clear-cutting; 
A2, 2 years after clear-cutting; A5, 5 years after clear-cutting; A8, 8 years after clear-cutting; A11, 11 years after 
clear-cutting. Bars mean standard errors. 


3.5 Age of recruits 


In the plots of 1 year before clear-cutting, the average ages of recruits were 8.44 (+4.49) and 7.38 
(+2.88) for EB1 and NB1, respectively. In the plots of 2 years after clear-cutting, the average ages 
were decreased to 5.86 (+1.08) and 7.05 (+1.74) for EA2 and NA2, respectively (Р<0.01 and 
P=0.142; Fig. 6). Then, the average age was gradually increased with time (Fig. 6). The age was 
influenced by the slope aspects. Before clear-cutting, the average age of recruits in EB1 plot was 
higher than that in NB1 plot (P=0.017; Fig. 6). However, after clear-cutting, the average ages of 
recruits in plots on north-facing slopes were higher than those on east-facing slopes (P<0.01, 
except for P=0.13 between ЕА11 and NA11; Fig. 6). In the plots of 2 years after clear-cutting, 
recruits originated from advanced regeneration accounted for 100% and 98.32% of total recruits 
for EA2 and NA2, respectively (Table 4). Then, the proportion of recruits originated from 
advanced regeneration was gradually declined with time (Table 4). The proportion of recruits was 
reduced to 33.33% and 38.05% for EA11 and NA11 plots, respectively (Table 4). 
Г] North-facing slope 
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Fig. 6 Variation of age of recruits in different slopes and plots. ВІ, 1 year before clear-cutting; A2, 2 years after 
clear-cutting; A5, 5 years after clear-cutting; A8, 8 years after clear-cutting; A11, 11 years after clear-cutting. Bars 
mean standard errors. 
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3.6 Spatial point pattern 

3.6.1 Univariate distribution patterns 

The spatial pattern of recruitment would be greatly influenced by time and slope aspect. The 
corresponding relations were analyzed with univariate analysis. The analysis revealed recruits 
were aggregated at small scales during most developmental stages (Fig. 7). On east-facing slopes, 
recruits were aggregated on the small scales (Fig. 7). For recruits of EB1, EA2, EAS, EAS, and 
EA11 plots, clumped pattern was observed on scales up to 9, 9, 6, 12, and 9 m, respectively (Figs. 
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Table 4 Characteristics of recruits originated from advanced regeneration after clear-cutting in different slopes 
and plots 


2 years after clear-cutting 5 years after clear-cutting 8 years after clear-cutting 11 years after clear-cutting 
EA2 NA2 EAS NAS EA8 NA8 ЕА11 МА11 
Quantitative 100.00 98.32 73.37 93.64 45.65 68.70 33.33 38.05 


proportion (%) 
Average height (m) 0.82+0.31 0.93+0.38 1.30+0.42 1.34+0.42 2.12+0.48 1.95+0.89 3.25+1.13 2.99+1.07 


Average age (а) 5.88x1.08 I MBE Loy) 7.11+1.36 8.16+1.67 10.83+1.59 10.97+2.02 14.70+2.63 14.05+1.96 


Note: Quantitative proportion was quantitative proportion of total recruits. Mean+SD. 


East-facing slope North-facing slope  Chronosequence 
1.67 (a) EBI 1.0 [ (b) NBI 
X 1 year before 
о clear-cutting 
0 5 10 15 20 25 30 
22 f (c) EA2 
~ 2.0 
> 2 years after 
o 18 clear-cutting 
1.6 
——S— 
Пааа ; 
0 5 10 15 20 25 30 0 5 10 15 20 25 30 
1.2 
= L0 5 years after 
o clear-cutting 
0.8 
T N 
Об 0.8 1 П n 1 1 ' 
0 5 10 15 20 25 30 0 5 10 15 20 25 3 
0.87 (g) EAS 
> 8 years after 
o lear-cutti 
5 clear-cutting 
0.4 1 1 1 1 1 П су уе ae eee 
0 5 10 15 20 25 30 0 5 10 15 20 25 30 
0.97 (i) EAM 
= 11 years after 
o { clear-cutting 
0.4 1 1 1 1 П J 3 
0 5 10 15 20 25 30 0 5 10 15 20 25 30 


Distance (m) Distance (m) 


Fig. 7 Univariate distribution patterns of recruits. The null model of complete spatial randomness (CSR) was 
applied in all 10 plots. The black lines indicate ring statistics O(r). The gray lines indicate the confidence 
envelope constructed by using the fifth smallest and largest values of 99 Monte Carlo simulations of the null 
model. O(r) that fall above, within, and below confidence intervals indicated clumped, random, and regular 
spatial patterns, respectively. 


7a, c, e, g and i). On north-facing slopes, spatial pattern of recruits was different during different 
developmental stages. For recruits of NB1 and МА? plots, clumped pattern was observed on 
scales up to 11 and 7 m (Figs. 7b and d). For recruits of NA5 plot, clumped pattern was observed 
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on very small scales (0—2 m), while recruits showed random pattern on other scales (Fig. 7f). For 
recruits of NA8 plot, random pattern was observed on scales 2—14 m while regular pattern on 
larger scales (Fig. 7h). For recruits in NA11 plot, the clumped pattern could be observed on all 
scales (Fig. 7j). On both slope aspects, the aggregation intensity was gradually decreased from 2 
to 8 years after clear-cutting (Figs. 7c—h). 

3.6.2 Bivariate association patterns 

The spatial pattern of recruitment would be simultaneously influenced by multiple factors. 
Bivariate analysis was applied to analyze the association of spatial pattern among different stages 
(recruits and canopy woods), as well as the association between seedlings and saplings. In EB1 
and NBI plots, recruits were spatially independent of canopy woods (Figs. 8a and b). 
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Fig.8  Bivariate association patterns between recruits and canopy woods. A bivariate antecedent condition null 
model was used in the 2 plots. The black lines indicate ring statistics O12(r). The gray lines indicate the 
confidence envelope constructed by using the fifth smallest and largest values of 99 Monte Carlo simulations of 
the null model. O12(7) that fall above, within, and below confidence intervals indicated attraction, independence, 
and repulsion, respectively. The description of gray and black lines is the same as Figure 9. 


Overall, spatial correlation between seedlings and saplings was transformed from attractive 
into independent over time. For the seedlings and saplings in EB1 and NBI plots, spatial 
attraction was observed on scales 0—6 and 0—11 m, spatial independence on scales 7—9 and 12-14 
m, respectively, and then it was shifted to spatial repulsions on all other scales (Figs. 9a and b). 
For the seedlings and saplings in EA2 plot, a spatial attraction was observed on scales up to 24 m 
and a spatial repulsion was observed on larger scales (Fig. 9c). However, for the seedlings and 
saplings in NA2 plot, spatial attraction was observed on scales 0—6 m, spatial independence on 
scales 7—13 m, and spatial repulsion on scales 14—30 m, respectively (Fig. 9d). For the EAS plot, 
seedlings were positively correlated to saplings on scales up to 28 m (Fig. 9e). In NA5, EA8, NA8, 
EA11, and NA11 plots, seedlings were spatially independent of saplings on all scales (Figs. 9f-j). 


4 Discussion 


4.1 Composition of tree species 


The natural recruitment was dominated by Chinese pine. The proportions of Chinese pine were 
higher than 90% in most of plots (Table 3). In addition, few seed sources of other native trees 
species (such as aspen, elm (Ulmus pumila), black locust, and Mongolian oak) could be observed 
nearby all stands. In Chinese pine plantations on the Loess Plateau of China, there was sufficient 
seed resource. However, the recruitments of other native tree species (including Liaodong oak 
(Quercus wutaishanica), Asian white birch (Betula platyphylla) and Amur maple (Acer ginnala)) 
would be limited by the thick litter covered on emerging seedling (Guo et al., 2015). Under the 
conditions of clear-cutting, the decomposition of litter would be accelerated and litter source 
would be reduced. The limitation from litter would be decreased over time. Then the seed sources 
would become a main limiting factor. This conclusion could be supported by other studies. And 
they found that if the clear-felled plantations were surrounded by the seed resource of other tree 
species, there would be relatively high proportion of other tree species (Yamagawa and Ito, 2006; 
Spracklen et al., 2013). 
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Fig. 9 Bivariate association patterns between seedlings and saplings. A bivariate independence null model was 
applied in all 10 plots. 


4.2 Density of recruits 


The density of recruits was significantly increased in the plots of 2 years after clear-cutting (Fig. 2). 
It was indicated that clear-cutting was conductive to seedling germination. Many environmental 
factors were correspondingly varied, including increased light intensity and soil temperature, as 
well as the declined thickness of litter layer and autotoxicity of Chinese pine (Li et al., 2010; Han 
et al., 2014). The gradual decrease of recruit density after clear-cutting was attributed to a high 
mortality of recruits (Fig. 2). The reason was the lack of protection from drought and heat after 
clear-cutting (Coates, 2002). What's more, severe competition for nutrients, water, and space 
would lead to self-thinning of recruits. 

Furthermore, slope aspect was observed to make remarkable effects on density of recruits. 
Before clear-cutting, the recruit density in EB1 plot was 56.59% higher than that in NBI plot (Fig. 
2). The tendency was consistent with previous results (Guo, 1992; Dong et al., 2013). For 
plantations with Chinese pine, the conditions were better on east-facing slopes than those on 
north-facing slopes, such as light and heat, which was more beneficial for germination and 
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survival of recruits (Guo, 1992). After clear-cutting, the water availability became a main limiting 
factor. In this area, the annual precipitation was limited (356.7—620.3 mm from 2004 to 2015; 
Table 1) and soil moisture content would be reduced after clear-cutting. After clear-cutting, the 
habitat was more suitable for recruit germination and survival on the north-facing slopes. 
However, the density of recruits was always higher on east-facing slopes after clear-cutting, 
except for EA5 and NAS plots (Fig. 2). This result was contrary to expectation. Thus, the reason 
of this phenomenon needed to be further explored. 


4.3 Height and DBH of recruits 


Both the quantity and proportion of saplings (higher than 0.5 m) in the plots of 2 years after 
clear-cutting were significantly higher than that in the plots before clear-cutting (Fig. 2; Figs. 
4a-d). In previous study, it was also found that the recruit height of Chinese pine could be 
promoted after clear-cutting (Liu et al., 2009). The development process of Chinese pine recruits 
can be divided into temperature-required phase and light-required phase. The seedlings of about 4 
years could be classified into light-required phase, thus the growth of saplings could be promoted 
with increasing light after clear-cutting (Liu et al., 2014). It was reported that the density, the 
average height and the average DBH of 11 years old Chinese pine plantations were about 
1500—2000 stems/hm?, 2.3—3.0 m and 2.1—2.8 cm on north-facing and northeast-facing slopes, 
respectively (Wang et al., 2012). In our study, the recruit densities (with the height 22.3 m) were 
1536 and 1725 stems/hm’, and the corresponding average DBH were 3.6 and 3.8 cm in EA11 and 
МА 11 plots, respectively. Therefore, we concluded that the speed of forest restoration was higher 
with naturel recruitment than that with artificial afforestation. 

In our study, recruit density and average DBH were not significantly varied in plots of 8 and 11 
years after clear-cutting (Figs. 2 and 5). Thus, the average DBH could not increase 
correspondingly with the increase of average height. In this interval, the self-thinning did not 
make obvious effects and the existing density had hindered the growth of saplings. Therefore, the 
artificial thinning was necessary in the plot of 8 years after clear-cutting. The proportion of 
recruits with higher height was dramatically declined in EA2 and NA2 plots compared to the EBI 
and NBI plots (Fig. 4). The results indicated a damage of advanced regeneration by the 
occurrence of clear-cutting. 


4.4 Age of recruits 


In our study, slope aspect had obvious effects on age of recruits (Fig. 6), and this result was 
consistent with the findings of previous studies (Han et al., 2012; Liu et al., 2014). In the plots of 
1 year before clear-cutting, the average age of recruits in NB1 plot was lower than that in EB1 
plot (Fig. 6). The limitation effects of light were more serious on recruits within light-required 
phase in NBI plot (Liu et al., 2014). Then it led to a lower average age. After clear-cutting, the 
average ages of recruits on north-facing slopes were always higher than those on east-facing 
slopes (Fig. 6). The reason can be attributed to the harsher habits and higher mortalities of recruits 
on east-facing slopes. The proportion of recruits originated from advanced regeneration was 
gradually declined with time (Table 4). Specially, the proportions of recruits originated from 
advanced regeneration were only 33.33% and 38.05% in EA11 and NAII plots, respectively 
(Table 4). However, the average heights of these recruits in both plots respectively reached 3.25 
(+1.13) and 2.99 (+1.07) m (Table 4). The results suggested that these recruits played vital roles 
in creating shade environment and improving microclimate. Above analysis suggested that 
advanced regeneration always played significant roles after clear-cutting throughout the stages of 
forest recovery. 


4.5 Spatial point pattern 


The main univariate distribution type of recruits in this study was aggregated at small scales (Fig. 
7). The aggregated distribution of natural recruitment has been observed in many studies 
(Camarero et al., 2005; Fajardo et al., 2006; Barbeito et al., 2008; Ignacio et al., 2009; Franklin 
and Santos, 2011). Both the favorable patches distribution and wind-spread seed character of 
Chinese pine were related to the aggregated distribution. The dispersal efficiency of wind-spread 
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seed would be intensively affected by various factors such as wind and surrounding trees. Thus, it 
would be conductive to forming aggregated distribution (Greene and Johnson, 1996; Condit et al., 
2000). The formation of aggregated distribution was helpful for forming close to natural forest. 
About 209%-50% of the forest land should be formed with open space for maximizing structural 
diversity and enhancing wildlife value (Rodwell and Patterson, 1994). The aggregation scales of 
recruits were no more than 12 m in all plots except for NA11 plot (Fig. 7). The seeds of Chinese 
pine were unable to be dispersed over a long distance, due to the lack of seed balloon, wings, and 
other auxiliary flying devices. Thus, the aggregation scale of recruits was limited. In the plots of 2 
to 8 years after clear-cutting, the aggregation intensity was gradually decreased because of the 
self-thinning effect (Figs. 7c—h). With the increased size of recruits, the competition for resources 
would also be increased among neighboring recruits. 

Recruits and canopy woods were independent of each other (Fig. 8). And other studies have 
found that before clear-cutting the environmental factors played more vital roles than 
density-dependent mortality (Franklin and Santos, 2011) and canopy wood nurse effects (Fajardo 
et al., 2006). When the stand age was low, recruits may be positively related to canopy woods. 
The reason may be the dispersal limitation and canopy wood nurse effect. However, with the 
increased competition between recruits and canopy woods, the positive correlation could be 
changed over time (Holmgren et al., 1997). The current spatial relationship would be finally 
formed. After clear-cutting, seedlings were positively correlated with saplings in a short time 
(Figs. 9c-e), indicating that the shade environment provided by saplings was beneficial to 
seedlings. 


5 Conclusions 


The density of advanced recruitment of Chinese pine plantations was extremely high with great 
potential. After clear-cutting, the potential of recruitment was quickly released. The recruits 
showed better growth statuses on north-facing slopes than on east-facing slopes after clear-cutting. 
Considering the density, height, and DBH of recruits, the speed of restoration with natural 
recruitment within 11 years after clear-cutting was faster than that with artificial afforestation. 
The height distribution of recruits indicated that the vertical structure diversity of naturally 
regenerated forest was higher than that of single-storied even aged plantation. The recruits 
established predominantly in clumps suggested that the horizontal structure diversity of naturally 
regenerated forest was also higher than that of plantation which was planted with uniform spacing. 
Thus, natural recruitment had the advantages of faster restoration speed and higher structure 
diversity than artificial afforestation in the early stage of forest recovery. The results indicated 
that natural recruitment may be possible to produce forest of a high economic, ecological and 
landscape value at a substantially reduced cost. Therefore, we proposed that clear-cutting of 
Chinese pine plantations recovered by natural recruitment would be a promising and effective 
method for establishing multifunctional forest. 

However, forest recovery was a long-time process while clear-cutting was only applied since 
2003, the observation should be continued to explore the temporal and spatial variation tendencies 
of natural recruitment. Several management strategies can be proposed in the light of our results: 
(1) artificial thinning was necessary in the plot of 8 years after clear-cutting; (2) advanced 
regeneration should be protected during clear-cutting operation, given its important role in early 
stage of forest recovery; and (3) artificial disturbance could be conducted to increase species 
diversity, such as implementing supplementary planting of other native tree species and retaining 
accompanying tree species in artificial thinning as far as possible. 
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